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Abstract 

The QCD analytic running coupling aan which has no nonphysical singularities 
for all > 0 is considered for the initial perturbation theory approximations up 
to four loop order. The finiteness of the analytic coupling at zero is shown to 
be a consequence of the asymptotic freedom property of the initial theory. The 
nonperturbative contributions to the analytic coupling are extracted explicitly. For 
all Q > A they are represented in the form of an expansion in inverse powers of 
Euclidean momentum squared. The effective method for a precise calculation of 
the analytic running coupling is developed on the basis of the stated expansion. 
The energy scale evolution of the analytic running coupling for the one- to four- 
loop cases is studied and the higher loop stability and low dependence on the quark 
threshold matching conditions in comparison with the perturbative running coupling 
were found. Normalizing the analytic running coupling at the scale of the rest 
mass of the Z boson with the world average value of the strong coupling constant, 
aaniM"^) = 0.1181 ± 0.002, one obtains as a result of the energy scale evolution of 
the analytic running coupling aan{M^) = 0.2943i'iQgjQg that is notably lower than 
the estimations of the coupling strength available at the scale of the mass of the r 
lepton. 
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I. INTRODUCTION 


The strong coupling constant Og is the basic parameter of Quantum Chromodynamics 
(QCD) and its determination appears to be one of the most important problems 
The perturbation theory supplemented with the renormalization group method works 
effectively beyond the infrared region. The nonphysical singularities of the perturbation 
theory arise in the infrared region of QCD and should be canceled by the nonperturbative 
contributions. The nonperturbative contributions arise quite naturally in an analytic 
approach to QCD (for a review see, e.g. [Q). The idea of the procedure goes back to 
Refs. [|, H devoted to the ghost pole problem in QFT. The foundation of the procedure is 
the principle of summation of imaginary parts of the perturbation theory terms. Then, the 
Kallen - Lehmann spectral representation results in the expressions without nonphysical 
singularities. In recent papers [0, H it is suggested to solve the ghost pole problem in QCD 
demanding the as(Q^) be analytic in (to compare with the dispersive approach 0). 
As a result, instead of the one-loop expression a^p{Q‘^) = (dvr/feo)/ ln(Q^/A^) taking into 
account the leading logarithms and having the ghost pole at = A? {Q^ is the Euclidean 
momentum squared), one obtains the expression 


a 


( 1 ) 


(Q^) = 


dvr 


+ 


A2 


ln(QVA2) A2-Q2 


( 1 ) 


Eq. (|l|) is an analytic function in the complex Q^-plane with a cut along the negative real 
semiaxis. The pole of the perturbative running coupling at = h? is canceled by the 
nonperturbative contribution [A^ ~ /i2exp{—47r/(6oaan(h^))} ^-t aan(h^) —^ 0] and the 
value ai^^(O) = dvr/feo appeared hnite and independent of A. The important feature of the 
” analyticization procedure” discovered SB is the stability property of the value of the 
’’analytically improved” running coupling constant at zero with respect to higher order 
corrections, Q!^^^(0) = Q!^^^(0) = a^^^(O). Though the derivative of the analytic running 
coupling is inhnite at zero, aan(Q^) turns out to be stable with respect to higher order 
corrections in the infrared region as a whole. 

The 1-loop order nonperturbative contribution in Eq. ( 0 ) can be presented as conver¬ 
gent at > A2 of constant signs series in the inverse powers of the momentum squared. 
For a ’’standard” as well as for iterative 2-loop perturbative input the nonperturbative 
contributions in analytic running coupling are calculated explicitly in Ref. |^. In the 


ultraviolet region the nonperturbative contributions can also be represented as a series 
in inverse powers of the momentum squared with different coefficients of the expansion. 
The nonperturbative contributions to aan(Q^) up to 3-loop order in analytic approach 
to QCD are studied in Refs. [O, |T^. To handle the singularities originating from the 


perturbative input the method which is more general than that of Ref. lOT was developed. 
In Ref. |T^ the momentum dependence of aan and its perturbative and nonperturbative 
components in the infrared region are analysed. For the standard perturbative input the 
higher loop stability and low sensitivity with respect to the c quark threshold matching 
conditions were found for aan- lu Ref. the nonperturbative contributions for the 


4-loop case are considered briefly. 

In this paper the momentum dependence of the analytic running coupling up to 4- 
loop order is studied. In parallel, the behavior of the perturbative component is given 
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for convenience of comparison. In Section 2 we generalize slightly the standard fonr-loop 
solution for hnd and study the spectral density for and then we prove the impor¬ 
tant property aa„(0) = dvr/fco- In Section 3 we extract from the initial perturbative 
contribution and hnd in an explicit form the nonperturbative contribution We 

develop a technique of integration in the vicinity of severe singularities of the perturbation 
theory in the infrared region and represent in the form of a hnite limits integral. In 
Section 4 the power series representation for at Q > A is obtained. In Section 5 we 
study the momentum behavior of aan- We consider it instructive to normalize aan and 
at Mz and then compare their behavior to estimate the nonperturbative contributions 
at all momenta. We consider two methods of matching of the solutions with different 
numbers Uf of active quark havors. Finally in Section 6 we give our conclusions. In the 
Appendix we give the explicit formulas which allow one to simplify the integration in the 
vicinity of the singularities of the standard perturbation theory input. 


II. FROM THE RUNNING COUPLING TO THE 
ANALYTIC RUNNING COUPLING 


The behavior of the QCD running coupling agiQ"^) is dehned by the renormalization group 
equation 




dasm 

dQ^ 


P{ois) — PoOil + 


( 2 ) 
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Here uj is the number of active quark flavors and C, is the Riemann zeta-function, 
C 3 = C(3) = 1.202056903.... The hrst two coefficients f3o, /3i do not depend on the 
renormalization scheme choice. The next coefficients do depend on it. Calculated within 
the MA-scheme in an arbitrary covariant gauge for the gluon held they appeared to be 
independent of the gauge parameter choice. Values of the coefficients (H) are given in 
Table |^. All these coefficients are small enough and decrease in absolute value with nj 
increasing. All the coefficients are negative except /32 at n/ = 6. 

The integration of Eq. yields 

^ ^ (/?o/^2 — {pi — 2j3oj3ij32 + PlPPj 

0(al(Q^)) = -/3oln(QpA^)+C. (4) 
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Table 1: Uf dependence of MS values of Pi {i = 0, 3), b, n, k. 



Po 

Pi 

P2 

P3 

b 

K 


0 

-0.87535 

-0.64592 

-0.71986 

-1.17269 

0.84298 

0.51033 

-1.16716 

1 

-0.82230 

-0.56571 

-0.58199 

-0.91043 

0.83663 

0.49541 

-1.20019 

2 

-0.76925 

-0.48550 

-0.45019 

-0.68103 

0.82045 

0.46922 

-1.26081 

3 

-0.71620 

-0.40528 

-0.32445 

-0.48484 

0.79012 

0.41467 

-1.36791 

4 

-0.66315 

-0.32507 

-0.20477 

-0.32222 

0.73920 

0.28506 

-1.56255 

5 

-0.61009 

-0.24486 

-0.09116 

-0.19354 

0.65784 

-0.07234 

-1.95343 

6 

-0.55704 

-0.16465 

0.01638 

-0.09914 

0.53061 

-1.33654 

-2.94623 


The integration constant is represented here as a combination of two constants A and 
C. Dimensional constant A is a parameter which dehnes the scale of Q and is used 
for developing the iteration procedure. Iteratively solving Eq. (|[) for at L = 

ln(QVA2) —>• cxo we obtain 




PoL + ^{\nL + C) 
Po 




2PIL^ 


{In L + Cf - 


PIL ' 
‘^PoP2 

~w 


(lnL + C')-l + 


PlP^ 


/Sf 


+ 0 


1 

1 ? 


(5) 


where C = ln(—/do) + {Pq/P i)C. Inverting Eq. (|]) one obtains 




PoL 




;3§i= 


^PoP2 


{InL + Cf - {InL + C) - 1 + 


(InL + Cf-plnL + Cf-h- 


n r I I ^ ^o/^3 



Within the conventional definition of A as Aj^ 


+ 0 

\L‘*. 

( 6 ) 

one chooses C = 0. At that the 
functional form of the approximate solution for turns out to be somewhat simpler, 

but it requires distinct A-j^g for different rif. With this choice, Eq. (|]) at the three loop 
level corresponds to the standard solution written in the form of the expansion in inverse 
powers of logarithms [|l| , and at the four loop level it corresponds to . We shall deal 
with nonzero C since this freedom can be useful for an optimization of the finite order 
perturbation calculations. Moreover, in the presence of the nj-dependent constant C it is 
possible to construct matched solution of Eq. (^) with universal nf independent constant 
A 0. 

Let us introduce the function a{x) = (6o/47r)as((5^), where x = Q’^ jIM. Then instead 
of (|) one can write 


a uc = 


Inx 


_ \n{\nx) + C 2 
In^x 


(ln(lna:) + C)^ ln(lna;) + C k 


Im X 


In'^ X 


+ 


In"^ X 


- ¥ 


(ln(lna;) + C)^ 5 (ln(lna;) + C)^ ln(lna;)+C' 

In"^ X 2 In^ a; In"^ a; 


K 


Im X 


• (7) 
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where the coefficients 6, k, and k are equal to 
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K 
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( 8 ) 


The values of parameters 6, k, and k of Eq. (|^) for different n/ are given in Table |l|. At 
X ~ 1 the perturbative running coupling is singular. At large x the 1-loop term of Eq. 0 
dehnes the ultraviolet behavior of a{x). However, for small x the behavior of the running 
coupling depends on the approximation we adopt and at x = 1 there are singularities of 
a different analytical structure. Namely, at x ~ 1 the leading singularities are 


a(^^(x) ~ —a^^^(x) ~^ln(x-l), 

X — 1 (x — Ij"' 

a^^\x) ~ ^ In^(x-l), a^^^(x)-- ^ ln^(x - 1). (9) 

(x — 1)-^ (x — 1)^ 


From Eqs. (^) we know the leading behavior at x ~ 1 of the additional terms which should 
cancel the perturbative singularities. But in principle it gives no information on their 
behavior at large x. The analytic approach removes all these nonphysical singularities in 
a regular way. 

The analytic running coupling is obtained by the integral representation 



( 10 ) 


where the spectral density p(cr) = Imaa„(—a — ffi). According to the analytic approach 
to QCD we adopt that Imaa„(—a — fO) = Ima(—a — iO), where a(x) is the perturbative 
running coupling. It is clear that dispersively-modihed coupling of form (llOl) has analytical 
structure which is consistent with causality. 

By making the analytic continuation of Eq. (0) into the Minkowski space x = —cr —iO, 
one obtains 

[ln(ln(T — in) + C]"^ 

(In cr — in)^ 


a(—a — iO) = 


In a 


-b 


In (In a — in) + C 2 


in 


(Incr 




in 


ln(ln(T — in) -|- C 


+ 


K 


(In a — in)^ (In a — in)'^ 


[ln(ln a — in) -|- C]^ 5 [ln(ln a — in) + 

2 (Incr — in)^ 


, , ln(ln(T — in) + C 

+ (3«: + 1)A^-+ 


(Incr — in)^ 
R 


(Incr —ivr)^ (Incr —f7r)^j 
Taking an imaginary part of Eq. o we hnd the spectral density 

p(cr) = p^^^(cr) -h Ap^^^(cr) Ap^^^(cr) Ap^^^(cr). 


( 11 ) 


( 12 ) 
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Ap^^^(cr) = 
62 


(62 + 7r2)2 
TT (36^ - vr2) (f2(6) - K 


27itFi{t) - {t^ - 7r2^ ^ 2 ( 6 ) 

— 26 (^62 — 37r2^ Fi(6)F2(6) 


(13) 

(14) 


(62 + 7r2)3 - 

— TT ^36^ — 7r2^ -^ 1 ( 6 ) +6 [t^ — F2{t) + ttk (^36^ — 7r2^ 


(15) 


Ap' 




^3 


^ 6 "^ - 67r262 + tt'^) (-Ff ( 6 ) - 3^2( 6 )^ 2 ( 6 )) + 47r6 {t^ - (^Ff{t) 

-Ff (6)) + 5 (^6"^ - 67r262 + tt'^) Fi{t)F2{t) 


(62 + 7r2)4 L 

-3Fi(6)F|(6)) - 107r6 (62 - 7r2) (^^2(6) 

+ 47r (1 + 3k) 6 (t^ — Fi{t) — (1 + 3^:) (6"^ — F2{t) + Ai^Rt (t^ — 

Here 6 = In (a), 


Fi( 6 ) = - ln( 6 ^ + TT^) + C, - 6 ^ 2 ( 6 ) = arccos 


-\/62 + 7r2 ’ 


• (16) 


(17) 


p*^^^(cr) is the 1-loop spectral density and Ap6)(cr) are higher loop corrections to the 
spectral density. With Eqs. ([I0|) , (|^) — (1T7|) the analytic running coupling can be 
studied, e.g. by numerical methods. For the 1 — 4-loop cases the spectral density of the 
analytic running coupling is shown in Fig. |I[ For the curves in Fig. ^ and in the next Fig. 
the parameter values C* = 0, nj = 3 are chosen. Beyond the 1-loop approximation one can 
see the higher loop stabilization of the spectral density and in the region of |6| > 10 it is 
practically the same for the 2 — 4-loop cases. Integrating the spectral density numerically 
with replacement of the inhnite limits by hnite cut parameter T leads to the relative error 
~ 1/T, and at large T it is important not to lose the higher loop contributions. In Fig. 
the higher loop corrections to the spectral density are shown. In fact we deal with rapidly 
oscillating functions and one needs special methods for a precise integration. F.g., for 
the 4-loop case it is difficult to get a 2-percent accuracy for aan{M^) using the standard 
integration program DGAUSS. 

We shall obtain another more effective method for precise calculation of aan{Q‘^) which 
is not connected with the numerical integration. 

Function a{x) in Fq. (|^ is regular and real for real x > 1. Thus, to hnd the spectral 
density p{a) we can use Schwarz reflection principle (a(x))* = a(x*) where x is considered 
as a complex variable. Then 

' -a — iO) — a{—a + 60)). 


p{(^) = Yi 

Let us introduce function 4*( 2 ;) of the form 


(18) 


<I)(2;) =- h -h 6 


(ln( 2 ;) -|- C) ln(z) + C ^ k 


- ¥ 


(ln(x)+C')^ 5(ln(^)+C')^ In(^) + G k 

-?- 2 -?-+ (3'‘ + l)—^ + 7 


(19) 
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To choose the main branch of the multivalued function ( p!9D we cut the complex z-plane 
along the negative semiaxis. Then solution (|^) can be written as a(x) = $(lnx). Function 
a(x) is unambiguously dehned in the complex x-plain with two cuts along the real axis, 
physical cut from minus inhnity to zero and nonphysical one from zero to unity. Then 


P(^) 


1 

Yi 


($(ln a 


ivr) — <F(lncr + m)). 


( 20 ) 


By the change of variable of the form a = exp(t), the analytical expression is derived 
from (0), (0 as follows: 


®an (^) 



dt -7 

X + 


X {<F(f 


m) — <F(f + ivr)} . 


( 21 ) 


Let us prove that aan(O) = 1. It follows from Eq. pl| ) that 


®an(0) 



dt — ivr) — <F(f + ivr)} 



— OO 


<F(f — ivr) 


1 - 
t — in. 


<F(f + ivr) 


1 ■ 
t + in. 


+ 


1 


t — in 


1 ■ 
t + in. 


( 22 ) 


For the hrst term in Eq. (p^ we close the integration contour in the lower half-plane of 
the complex variable t by the arch of the ’’inhnite” radius without affecting the value of 
the integral. We can do it because the integrand multiplied by t goes to zero at | f |—>• cxd. 
There are no singularities inside the contour, and thus we obtain a zero contribution from 
the term considered. For the second term we close the integration contour in the upper 
half-plane of the complex variable t with the same result. Therefore we have: 


®an(0) 



— OO 


■ 1 
.t — in 


1 ■ 
t + in. 


1 . 


(23) 


For any hnite loop order the expansion structure of the perturbative solution in inverse 
powers of logarithms ensure the property of the analytic coupling aan(O) = 1. The argu¬ 
ments are suitable for all solutions as long as the singularities are situated at the 
real axis of the complex 2 ;-plane, in particular for the iterative solutions of Refs. 0, |^. 


III. EXTRACTION OF THE NONPERTURBATIVE 
CONTRIBUTIONS 

Let us see what the singularities of the integrand of (pT]) in the complex t-plane are. First 
of all the integrand has simple poles at f = Inx ± in{l + 2n), n = 0, 1,2,.... All the 
residues of function exp(t)/(a; -|- exp(f)) at these points are equal to unity. Apart from 
these poles the integrand of (pl|) has at f = ±in poles up to fourth order and logarithmic 
type branch points which coincide with the poles from the second order to the fourth 
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order. The initial integration contour in the complex t-plane and singularities of the 
integrand of Eq. (|2l|) are shown in Fig. 3(a). Let us cut the complex t-plane in a standard 
way, t = iivr — A, with A being the real parameter varying from 0 to oo. Further on we 
append the integration by the arch of the ’’inhnite” radius without affecting the value of 
the integral, and close the integration contour Ci in the upper half-plane of the complex 
variable t excluding the singularities at t = in. In this case an additional contribution 
emerges due to the integration along the sides of the cut and around the singularities at 
t = in. The corresponding contour we denote as C 2 . The integration contour Ci is shown 
in Fig. 3(b), whereas the contour C 2 is shown in Fig. 3(c). 

Let us turn to the integration along the contour Ci. For the integrand of Eq. (0. 
which we denote as F(t), the residues at t = lnx-|-i7r(l -|-2n), n = 0,1, 2,... are as follows: 

ResF(f) \t=inx+iwii+ 2 n)= <h(lnx 2nin) - $(lna; 2ni{n + 1)). (24) 

By using the residue theorem one readily obtains the contribution S(x) to the integral (^I|) 
from the integration along the contour Ci 

1 r °° 

S(a;) =- / F{t) dt = 'Y] ResF (t = Inx -|- in{l + 2n)) = <h(lna;) = a{x). (25) 

2ni J 


Cl 


One can see that this contribution is exactly equal to the initial Eq. (|^). Therefore we 
call it a perturbative part of aan{x), a^^{x) = S(x). The remaining contribution of the 
integral along the contour C 2 can naturally be called a nonperturbative part of aan{x), 


aan{x) = a’“{x) + (x). 


(26) 


Let us turn to the calculation of a^^{x). We can omit the terms of the integrand in 
Eq. (0) which have no singularities at t = in. Then we have 


(^) = 

C 2 
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X -\- \t — in 


-h 


ln(t — in) + C 


{t — in)"^ 


- F 


+F 


[ln(t — in) -1- C]'^ ln(t — in) + C 


{t — in)^ 
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[ln(t — in) + CY 5 [ln(t — in) -\- C\ 


{t — in)‘^ 


{t — in)* 


-|- (3k 1) 


ln(t — in) + C 


+ 


K 


(t — in)* (t — in)* 


Let us change the variable t = z + in and introduce the function 

1 


fiz) = 


1 — a; exp (— 2 ;) 


(27) 

(28) 


Then we can rewrite Eq. (|27|) in the form 
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Ox) = 


2ni 


[ dzf(z) | - - 6 

rin(z) 01 

+ F 

J [z 
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2 ;^ 2 ;^ 



In^(z) 


+ (2C-1)!H + 
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c + c^ 


7 


























5\ In' 


+ SC' - - — + (3C" - 5C + 3k + 1)— + 


Inz C3 - ^ (3k + 1)C + K 


ln3(z) 


(29) 

The cut in the complex z-plane goes now from zero to —oo. Starting from z = —cxo — iO, 
the contour C goes along the lower side of the cut, then around the origin and then 
further along the upper side of the cut to z = — oo + iO. a: is considered here as a real 
variable, x > 1. Then the contour C can be chosen in such a way that it does not 
envelop ’’superfluous” singularities, and the conditions used in the Appendix for hnding 
the corresponding integrals are satished. Function (pSD with its derivatives 


xexp(-z) xexp(-z)(l+ xexp(- 2 ;)) 


nz) = 


(1 — xexp(— 
xexp(— 2 ;) 

(1 — xexp(—; 2 ))^ 


(1 — xexp(— 

+ 4xexp(—z) + exp(— 22 ;)^ , 


(30) 


/'"'(z) = A ^ llxexp(-; 2 ) + llx^exp(-2^) +x^exp(-3z)) 

(1 - xexp(-;2))5 V ^ vv y vv jj 

decrease exponentially at z —*• — 00 . Therefore, we shall omit the boundary terms in 
formulas given in the Appendix!. Then, from Eq. (^9]) one can obtain 


- —Jdz!^f{z)ln{z)-b 


{l + C)ln{z) + -ln\z) 


nz) 
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-b‘^ 

2 
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(2 + K + 2C + C^) \n{z) + (1 + C) \n^{z) + - \n%z) 


riz) 


- 6 ^ \(Q + —K + K + 3{2 + K)C + 3C^ + cAlnz + -{2 + K + 2C + C^)ln'^z 
6 LV 2 / 2 


+ {l + C)\n^z + -hi^z 
4 . 


/""(*)} . 


(31) 


Taking into account that function f{z) with its derivatives is regular at real negative 
semiaxis of 2 ; we can rewrite equation (|^ in the form 


«rw= - 
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1 . 


du \ f'{u)Ai{u) — b 


(1 + C)Ai(m) + -A2(u) 


/"(«) 


+ K + (1 + C)^^ Ai(u) + (1 + C)A2 ('u) + —A3(m) 


/"'(«) 


^2 + —K + K + 3(1 + C)(l + k) + (1 + C)^^ Ai('u) 


- 

- 

+ 2 ^ C)^^ A2(m) + (1 + C)A3(-u) + —A4(-u) 

where u is real, n < 0 and Ai{u) are discontinuities of the powers of the logarithms 


r{u) , (32) 


^i(^) = 7 ;—: (ln(M + iO) — ln(M — iO)) = 1, 
2tii 


^Using these formulas with / = 1 one can make sure that a{x = 0) = 1. The boundary terms should 
be considered in this case. 
























( 33 ) 


A 2 (m) = -+ iO) — — iO)) = 2 ln(— m), 

27ri ^ ' 

A 3 (m) = 77 —: (ln^(M + iO) — ln^(M — zO)) = 3ln^(— m) — TT^, 

27rz ^ ' 

A 4 (m) = (ln^(M + zO) — ln'^(M — zO)) = 41n^(—Zi) — 47 r^ln(—zz). 

we obtain 


Let us introduce the variable a = exp(zz). From Eqs 
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TT 
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X ln(— Incr) + 3(1 + C) ln^(— Incr) + ln^(— Incr) 
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^ + llx^cr + llxcr^ + cr^ 


(x — cr)5 

Integrating the terms of Eq. (|3^ independent of logarithms one can obtain 
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= -- 
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1 l(l + C')x 
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x(x + 1 ) 
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TT 
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x(x^ + 4x + 1) 
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1 2 

+ X J da 3^1 —^ + zz + (l + Cy^ ln(— In cr) 


+ 3(1+ C')ln2(-Incr) +ln3(-Incr) 


X'^ + llx^cr + llxcr^ + cr'^ 


X - cr 


(35) 


This formula gives the nonperturbative contributions in an explicit form. 


IV. BEHAVIOR OF THE NONPERTURBATIVE 
CONTRIBUTIONS AT Q > A 

Let us turn to the large Q behavior of the nonperturbative contributions. The following 
expansions appear to be useful (x > 1 > cr > 0 ) 


X — 1 


00 1 

= E 


X 


_ zz x(x + c^) _ 

r^n’’ t rr. 


00 ^ 2 ^n-l 
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(x - 1)3 ^ = Z- 
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[X - a) 
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x{x^ + 4x + 1) 
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oo ^3 
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E l L 

iy>7X ^ 


vx x(x^ + llx^a + llxcr^ + cr^) n’cr" 


cxD n—1 


n=l 


(x — aY 


= E 

n=l 




Note that the coefficients in Eqs. (|36|) are monomials in powers of n. Expanding Eq. 
in the inverse powers of x with using Eqs. (^) we have 

“S‘(^) = E (37) 


n=l 


X’‘ 


where 


Cn = —I + bn ll + C + n J da a"' ^ In (— hi(cT)) 


2(1 + C) In (— In(cr)) + In^ (— In(cr)) 


--6^77,^ <1 + K-l■(l + + n J da a'^ 

{ 0 

|2 + + /t + 3(1 + C) ^1 - Y + (1 + CY + n J da a"'~^ 3^ 

+ (1 + ln(— In a) + 3(1 + C) ln^(— In a) + ln^(— In a) |. 

Making the change of variable a = exp(—f) and integrating [^, over t one can End 

1 oo 

y dcr In (—In(cT)) = J dte~'^^\n{t) = -(ln(n) + 7 ), 


TT^ 

1-- + K 


(38) 


/'<i(Tcr”-Mn7-lii(cr)) = f (ii e"”'lii^(i) = 

0 0 
1 00 

J (icrcr"'“Mn^ (—In(cr)) = J dt e~"‘^ liY{t) = —— 


(ln(n) + 7 )^ + ^ 
0 


(39) 


(ln(n) + 7 )^ + Y(ln(n) + 7 ) + 2(3 


Here 7 is the Euler constant, 7 ~ 0.5772. From Eqs. (^), (^) we finally have 


2„2 


Cn = —1 + [1 + (7 — 7 — ln(? 7 ,)] — -^b~n 


l-^ + n+(l + C- ^- ln(n)) 


-b^n^ 


2 + + K - 2(3 + (1 + C - ^ - ln(n)) 


TT 


7 - ln(n)) [1- — + K 


(40) 


We can see from Eq. (^) that power series (|^) is uniformly convergent at x > 1 and 
its convergence radius is equal to unity. The resulting Eq. (|40D is scheme independent in 
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the sense that nj dependence is not hxed here, and the method used above allows one 
in principle to calculate next loops contributions to clarify the general structure of the 
coefficients c„. 

For numerical evaluation of the coefficients c„ we choose the MS scheme values of 
K, R and assume that (7 = 0. Then the coefficients c„ are dependent on n, Uf, and 
on the number of loops taken into account. In Table ^ we give the values of c„ and 
loop corrections for nj = 0,3,4, 5,6. The 1-loop order contributions to Cn are equal to 
— 1 for all n and nj. Up to 4-loop approximation the coefficients c„ for all n, Uf are 
negative. With the exception of the 3-loop case at n/ = 6, the 2 — 4-loop coefficients c„ 
for Uf = 0, 3,4, 5, 6 monotonously increase in the absolute value with increasing n. In the 
ultraviolet region [x ^ 1) the nonperturbative contributions are determined by the hrst 
term of the series One can see that for all n/ up to four loops ci is of the order of 

unity. The account for the higher loop corrections results in some compensation of the 
1-loop leading at large x term of the form 1/x. 


V. MOMENTUM DEPENDENCE OF 


The expansion coefficients increase in the absolute value not too fast and therefore the 
representation of the analytic running coupling of QCD in the form 


dvr °° 

aan{Q^) = aP\Q^) + — J2c, 

n=l 



(41) 


with Cn as in Eq. (^0|) provides one with the effective method for the calculation of aan at 
Q > A. At that there is no need for the summation of large number of terms of the series. 
Let us see what the convergence properties of the series (|37|) are. Since pn > In^(n) for 
all n > 1 and p > Pq = (3/e)^ ~ 1.4, one can consider the series 


54 


E 

n=l 





x{x^ + llx'^ + llx -I- 1) 
(x — 1)^ 


(42) 


as a comparison series for Eq. (|37| ) with coefficients (^) 
the series Eq. (^^ are not worse than that for the series 
A^-terms approximation of the series 


The convergence properties of 
^). The absolute error for the 


IS 


Af) = 


x^{x — 1 ) 


x{x^ -f llx^ -I- llx -I- 1) iNx^x"^ -I- 4a; -1- 1) 


{x — 1)^ 

6N^x(x + 1) AN^x 
-^-L 

X — 1 


(x — 1)^ 


+ 


-f W 


(x- 1)= 


(43) 


It is dependent on Q and N (with given nj, A) and the larger Q and N are, the smaller 
it is. For rather small x = 2 {Q = 1.4A) from Eq. (^3l) we hnd that the error of the 
approximation of the series (|4^ ) for A = 50 is ~ 10“® and for N = 100 it is ~ 10“^^. For 
larger Q there are no reasons to sum a large number of terms. For the approximation of the 
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Table 2: The dependence of Cn and loop corrections on n and Uf for the 1 — 4-loop cases. 




n 

^1—loop 

\ 2—loop 

A S —loop 

A A —loop 

^2—loop 

^3—loop 

^A—loop 


= 0 

1 

-1.0 

0.35640 

-0.01568 

-0.03900 

-0.64360 

-0.65929 

-0.69828 


2 

-1.0 

-0.45582 

0.08741 

-0.16455 

-1.45582 

-1.36841 

-1.53296 



3 

-1.0 

-1.70912 

-1.03012 

-0.89283 

-2.70912 

-3.73924 

-4.63207 



4 

-1.0 

-3.24886 

-4.51236 

-5.11707 

-4.24886 

-8.76122 

-13.87829 



5 

-1.0 

-5.00160 

-11.31238 

-18.56032 

-6.00160 

-17.31398 

-35.87430 



6 

-1.0 

-6.92407 

-22.24971 

-49.77660 

-7.92407 

-30.17378 

-79.95039 



8 

-1.0 

-11.17217 

-59.34790 

-213.31981 

-12.17217 

-71.52007 

-284.83987 



10 

-1.0 

-15.84626 

-120.76945 

-616.88776 

-16.84626 

-137.61571 

-754.50348 


= 3 

1 

-1.0 

0.33405 

0.01608 

-0.07825 

-0.66595 

-0.64987 

-0.72812 


2 

-1.0 

-0.42724 

0.19624 

-0.37379 

-1.42724 

-1.23101 

-1.60480 



3 

-1.0 

-1.60196 

-0.63626 

-1.28115 

-2.60196 

-3.23823 

-4.51937 



4 

-1.0 

-3.04517 

-3.48651 

-5.07338 

-4.04517 

-7.53168 

-12.60506 



5 

-1.0 

-4.68801 

-9.19185 

-16.30462 

-5.68801 

-14.87987 

-31.18449 



6 

-1.0 

-6.48996 

-18.47225 

-41.82403 

-7.48996 

-25.96221 

-67.78624 



8 

-1.0 

-10.47171 

-50.22832 

-174.16411 

-11.47171 

-61.70003 

-235.86414 



10 

-1.0 

-14.85275 

-103.11451 

-499.79465 

-15.85275 

-118.96725 

-618.76190 


= 4 

1 

-1.0 

0.31252 

0.04949 

-0.11006 

-0.68748 

-0.63799 

-0.74805 


2 

-1.0 

-0.39970 

0.31341 

-0.52880 

-1.39970 

-1.08630 

-1.61510 



3 

-1.0 

-1.49872 

-0.23818 

-1.51417 

-2.49872 

-2.73690 

-4.25107 



4 

-1.0 

-2.84891 

-2.48499 

-4.77483 

-3.84891 

-6.33389 

-11.10872 



5 

-1.0 

-4.38587 

-7.15989 

-13.83276 

-5.38587 

-12.54576 

-26.37852 



6 

-1.0 

-6.07168 

-14.89306 

-34.04904 

-7.07168 

-21.96474 

-56.01377 



8 

-1.0 

-9.79681 

-41.69613 

-138.28745 

-10.79681 

-52.49294 

-190.78040 



10 

-1.0 

-13.89549 

-86.71011 

-394.96378 

-14.89549 

-101.60560 

-496.56937 


= 5 

1 

-1.0 

0.27813 

0.11653 

-0.16002 

-0.72187 

-0.60535 

-0.76537 



2 

-1.0 

-0.35571 

0.55755 

-0.75021 

-1.35571 

-0.79817 

-1.54837 



3 

-1.0 

-1.33377 

0.50736 

-1.78434 

-2.33377 

-1.82641 

-3.61075 



4 

-1.0 

-2.53536 

-0.73077 

-4.12859 

-3.53536 

-4.26613 

-8.39472 



5 

-1.0 

-3.90317 

-3.73728 

-9.82126 

-4.90317 

-8.64046 

-18.46172 



6 

-1.0 

-5.40344 

-9.01126 

-22.11892 

-6.40344 

-15.41470 

-37.53362 



8 

-1.0 

-8.71859 

-28.07388 

-85.51994 

-9.71859 

-37.79247 

-123.31240 



10 

-1.0 

-12.36617 

-60.94080 

-243.69211 

-13.36617 

-74.30698 

-317.99908 


= 6 

1 

-1.0 

0.22433 

0.25378 

-0.22731 

-0.77567 

-0.52189 

-0.74920 



2 

-1.0 

-0.28692 

1.07460 

-1.01673 

-1.28692 

-0.21231 

-1.22904 



3 

-1.0 

-1.07581 

1.93179 

-2.00536 

-2.07581 

-0.14402 

-2.14938 



4 

-1.0 

-2.04500 

2.37205 

-2.96183 

-3.04500 

-0.67295 

-3.63479 



5 

-1.0 

-3.14826 

2.01775 

-4.07318 

-4.14826 

-2.13052 

-6.20370 



6 

-1.0 

-4.35837 

0.54419 

-6.02095 

-5.35837 

-4.81418 

-10.83512 



8 

-1.0 

-7.03234 

-6.87467 

-17.72700 

-8.03234 

-14.90701 

-32.63400 



10 

-1.0 

-9.97445 

-21.85073 

-53.77987 

-10.97445 

-32.82518 

-86.60506 
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analytic running coupling with only one first term of the series (|^) for the nonperturbative 
contributions taken into account, 


0.(1 


,{Q^) ~ a^W) - ^ |l - 6(1 - 7) + it' fl 


71 


- + K + (1 


7 )' 


— 

6 


2 + -K + K — 2(^3 + (1 — 7)^ + 3 (1 — 7) fl “ 


6 


+ K 


(44) 


the relative error was studied in Ref. |]^ for the 1 — 4-loop order cases. The approxi¬ 
mation of the nonperturbative ’’tail” by the leading term has been shown to give a one 
percent accuracy for Uan already at Q ~ 5A. 

In Fig. 4 the x dependencies of ttan, ctan presented for the 1 — 4-loop order 
cases. The nonperturbative contributions have been calculated by the series summation 
and the analytic running coupling has been calculated through the dispersive represen¬ 
tation (pRl). It turned out that the numerical integration in the cases considered expects 
dehnite caution. Insufficient accuracy of integration can look as an ungrounded stability of 
the analytic running coupling behavior with increase of the order of approximation. The 
equality (with the accuracy of 2 ■ 10“^ percent) of aa„(x) calculated through the dispersive 
representation and the sum of a^^{x) and calculated as the series for all x from 2 to 
20 served us as a criterion for the integration precision. The perturbative component 
increases with the decrease of x reaching unity at x 3(Q 1.7A). The nonperturbative 

component is negative (at a: > 1), it decreases with x compensating for the increase of the 
perturbative component. According to representation m the quantity aan{x) is regular 
for all a; > 0 and a£)^(0) = 1 (/ is the number of loops of the approximation). Though 
the derivative of aan{x) is infinite at zero we however make sure numerically of the higher 
loop stability of aan{x) in the infrared region. As seen in Fig. 4, the 3-loop and 4-loop 
analytic curves practically coincide even before the normalization at some finite point. As 
for the corresponding perturbative curves which have no common point at zero, they are 
not close to each other already at a: < 5. 

Let us consider the momentum dependence of aan (and for comparison) in the low 
momentum region provided that all solutions are normalized at the central point of the 
world average value a(M|) = 0.1181 ±0.002 |^. In this case the heavy quark thresholds 
should be taken into account. It seems natural to demand the analytic running coupling be 
continuous across thresholds. Let us adopt for aan, and for all 1 — 4-loop order cases 
the normalization condition = 0.1181, Mz = 91.1882 GeV and the matching 

conditions ruh = 4.3 GeV and 

rric = 1.3 GeV I. The corresponding sets of parameters A are given in Table 3. As seen 
in Table 3, ~ since the nonperturbative contributions in fact die out at 

the scale of normalization i. The momentum dependence of aan-, for the 1 — 4-loop 
order cases is presented in Fig. 5. As seen from Fig. 5, the 2 — 4-loop curves for aan 


•^The matching conditions sensitivity of aan will be considered slightly later. 

"^If to normalize the solutions as = 0.35, Mr = 1777.03 MeV § with the same matching 

conditions one obtains substantially larger values of Aan (e.g., for the 4-loop case ~ 630 MeV, 

490 MeV, 350 MeV, respectively). At that the higher loop stability is also observed, for the 2 — 4-loop 
cases ) = 0.128. 
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Table 3: The parameters (MeV), (MeV). Tif is the number of active quark 

flavors, the number of loops is indicated. The normalization and matching conditions are 
aK=5)(M|) = 0.1181, Mz = 91.1882 GeV, rrih = 4.3 GeV, 

GeV. 



1-loop 

2-loop 

3-loop 

4-loop 


=3) 

143.77 

372.50 

328.98 

332.50 

aV/= 

=4) 

120.55 

325.91 

289.67 

291.39 

aG/= 

=5) 

88.35 

227.51 

209.54 

209.53 

aG/= 

=3) 

150.64 

454.21 

382.30 

389.50 

aG/= 

=4) 

121.61 

339.90 

298.91 

301.64 

aV/= 

^ ^an 

=5) 

88.35 

227.60 

209.60 

209.61 


Table 4: The values with = 1777.03 MeV. The nor¬ 
malization and matching conditions are = 0.1181, Mz = 91.1882 GeV, 

Q-G/=5)(^/7i2^ = = 4.3 GeV. 



1-loop 

2-loop 

3-loop 

4-loop 


0.2740 

0.2930 

0.2943 

0.2943 


0.2802 

0.3262 

0.3179 

0.3230 


practically coincide indicating the higher loop stability of the analytic running coupling. 
The corresponding values of for the 1 — 4-loop cases are given 

in Table 4. 

Note that in Table 4 the values of 3-loop and 4-loop aan{M^) are the same and the 4- 
loop apt{M^) coincides with as{Mr) of Ref. 0 from r decays. Thus, the extrapolation to 
the energy scale Mz using the 4-loop solution for as with 3-loop matching at the bottom 
quark pole mass Mi, = 4.7 GeV made in Ref. results in as{Mz) = 0.1181. This value 
is used in our calculations. 

Let us consider the threshold matching conditions sensitivity of We fix (Af|) 

at its world average value and vary the value of the matching point corresponding to 
the 6-quark threshold. Then from the matching condition we 

can hnd the dependence of the parameters on the matching point value. For a 

rather wide interval of /i^ it is shown in Fig. 6. We see that parameters for the 

analytic running coupling go higher than that for the perturbative coupling. In the region 
of 4 - 5 GeV the /ift-dependence of A^an~^'^ is somewhat weaker than that of The 

dependencies of aan{M^), a^^^M^) on the matching point for the 1 — 4-loop order 
cases are shown in Fig. 7. For the analytic coupling the curves go lower than the corre¬ 
sponding curves for the perturbative coupling. The analytic coupling is much more stable 
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than the perturbative one with respect to higher loop corrections. In the region of 4 - 
5 GeV the /ib-dependence of aan{M^) is considerably weaker than that of In 

particular, for a{M‘^) = 0.1181, mi, = 4.3 ± 0.2 GeV 


aan{M^) = 0.29431 


+0.0004 — n 

0.0003) “ ~ 'J-'J^'J'J-O.OOOS- 


(45) 


We give here the results for the 4-loop aan and 

Let us consider one more heavy quark threshold matching condition. In the framework 
of the perturbation theory there is the prescription to connect the couplings with 
different nj according to which the coupling can be discontinuous at the matching point 
Hh- The idea of an implementation of this nontrivial matching conditions is to make the 
results (e.g., the connection between a{M^) and q;(M|)) be not substantially dependent 
on the exact value of the matching point |^. This conditions take the most simple form 


for two cases. First, = 'mh = mh{mh) where mh{fi) is the running MS mass of the 
heavy quark and second, pb = Mh with Mh being the heavy quark pole mass. Ghoosing 

m 


the hrst one according to 


we have 


for the 1-loop and 2-loop cases. 




1 + C2(a^ 


for the 3-loop case, and 


for the 4-loop case. Here 
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564731 82043 


72’ 124416 


27648 


Cs 


) 

(46) 


(47) 

y+C3(^a^^f\ni)/7iy 

(48) 

2633 , 

(49) 

31104*^’^-^ 


o 


correspond to Ref. 0. The 


The notations for the coefficients in Eqs. 
coefficients c„ of the present paper have its own dehnition. For this variant of matching 
conditions the dependencies of the parameters AG/=^) on the matching point value are 
close to those shown in Fig. 6 and we do not give the corresponding hgure (according 
to Eq. (^) both matching methods give the same results for the 1, 2-loop order cases). 
The dependencies of on the matching point /if, for the 3, 4-loop order 

cases are also close to the previous case of the continuous matching shown in Fig. 7. For 
a(M|) = 0.1181, TUf, = 4.3 ± 0.2 GeV 


aaniM^r) = 0.2947^™, = 0.32351°; 


0008 

0004- 


(50) 


The results are given for the 4-loop case. It is seen from Eqs. (|45D , (^) that the values 
of aan{M^) for two methods of matching are very close to each other (this is true also for 
aP\M^)). 

In Fig. 8 Q(an(Q^), q;^*(Q^) are shown for the 1 — 4-loop order cases with the normal¬ 
ization condition = 0.1181, Mz = 91.1882 GeV and continuous matching at 

mb = 4.3 GeV, me = 1.3 GeV. Without going into the details, we give in this hgure the 
data from Table 6 of Ref. [Q for the world summary of measurements of 
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VI. CONCLUSIONS 


In contrast to the recent papers we apply the analytic approach to the perturbative 
QCD running coupling constant in the form of the standard expansion in the inverse 
powers of logarithms up to the four loop order. An introduction of the complex variable 
t in the spectral representation for the analytic running coupling and study of the singu¬ 
larities structure of the integrand allowed one to divide the analytic running coupling into 
perturbative (initial) component and nonperturbative one (appeared as a consequence of 
’’forced” analyticity) exactly, as it is illustrated in Fig. 3. These components turned out 
to be connected with different singularities in the complex f-plane. It is shown that the 
nonperturbative contributions can be represented in the form of the expansion Eq. (^7[) in 
inverse powers of the momentum squared where the coefficients c„ are dehned by Eq. (^Ol). 
Eq. (1^) gives the effective method non-connected with numerical integration for calcula¬ 
tion of the analytic running coupling at Q > A with the calculation accuracy of standard 
mathematical functions. It can be important for making popular the considered variant 
of aan- In practice, for Q corresponding to n/ = 5, it is sufficient to take account of the 
leading nonperturbative term, as in Eq. (|4^) . 

On the basis of the developed method we study the momentum dependence of aan 
giving at the same time the behavior of the perturbative running coupling. To £x the 
solutions we used for all of them the same normalization condition at Mz where the 
nonperturbative contributions are negligible quantities. We can see in Fig. 4, Fig. 5, 
Fig. 7, and Fig. 8 the higher loop stability of the analytic running coupling for all Q > 0 
(the 1-loop case falls out of the common picture). For the perturbative case the higher 
loop stability takes place only at sufficiently large Q. 

We considered two variants of heavy quark threshold matching conditions for aan- 
The results appeared to be very similar. We showed the 6-quark threshold matching 
conditions stability of the analytic running coupling H. As a criterion we considered the 
dependence on the matching point fib of the correspondence of aan{M^) to aaniM^) for 
the 1 — 4-loop cases (for comparison a^^ was considered simultaneously). The situation 
is illustrated by Fig. 7. The energy scale evolution of the analytic running coupling gives 
a{M^) = 0.2943lo;oQ03 for fhe normalization at the world average value of a{M‘^) = 0.1181 
and matching by continuity with = 4.3 ± 0.2 GeV. With the same normalization 
condition at the scale of Mz for both matching methods considered aan{M^) is about 
0.03 less than a^^{M^). Therefore, if one regards aan as a true running coupling constant 
the noticeable discrepancy with r lepton decay data ii arises. 

A possible solution of this problem can be found if one changes the normalization 
condition at Mz- Let it corresponds to the value of Refs. [^, ^ a{M‘^) ~ 0.124 which 
is appreciably larger then the conventional one. Then for the 4-loop case the result is 
aaniM^) = 0.3270]'^o;ooo 3 with continuous matching at rub = 4.3 ± 0.2 GeV. 

As seen in Fig. 8, the analytic approach gives the running coupling which does not devi¬ 
ate essentially at sufficiently large momentum values from the usual perturbative running 
coupling constant. In the infrared region this approach allows one to solve the principal 
difficulty connected with nonphysical singularities. The question arises whether the ap¬ 
proach described takes into account the nonperturbative contributions to the right degree. 


^The c-quark threshold matching was considered in Ref. ha with same result. 
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There is a whole series of the approaches in which the nonperturbative contributions to 
the Green functions and running coupling are studied. These approaches are beyond the 
scope of the present paper and we only point out some papers |^, |^, |^, |^, |^, 
dealing with the nonperturbative contributions to the running coupling. 

To summarize, the analytic running coupling seems to be a good basis for the problem 
of’’genuine nonperturbative” contributions in ’’physical” as, which needs further analysis. 
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APPENDIX 


We give here the identities we need in our computations which can be obtained by means 
of an integration by parts. Let function f{z) of complex variable 2; be regular in some 
domain D where z = 0 E D. Dealing with the singularities of the integrands at the origin 
of the pole type coinciding with the logarithmic type branch points we cut the domain 
D along real negative semiaxis. Then for any contour C in the cut domain D which goes 
from 7^ 0 to ^2 7^ 0 one can End 

/ —f{z) = - f dz \n{z)f\z) + \n{z)f{z) 

J_ z J 21 

c c 


c c 

/ J/(^) = \ Hz)f{z) 


c 


c 


= -6/“^" ln(.-)/""(.-) + {-3pj/(.-) - - ^J"U) + gta(^)r(.-) 

c c 

J I \n(z)f(z) = -\jdz \n\z)f(z) + i \n\z)Hz) 


c 


c 


J J Hz)f{z) = - J dz (^ln(2;) + ^ ln^(2;)^ f'(z) + l~fiz) - ^-^f{z) + \n{z)f'{z) 
c c ^ 
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J z 
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c 


+ 


216^2 



575 ln(z) 
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Figure 1: The spectral density of the analytic running coupling up to four loop order. 














Figure 3: Complex t-plane integration. Perturbative contributions arise from the poles 
at t = Inx + f7r(l + 2n), n = 0,1, 2,.... Nonperturbative contributions emerge from the 
singularities at f = in. 
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Figure 4: The analytic running coupling Uan and its perturbative component and 
nonperturbative component as functions of a; = Q‘^/A‘^ for the 1 — 4-loop order 

cases. Here Uf = 3. 








Figure 5: The momentum dependence of aan, for the 1 — 4-loop order cases. The 
normalization and matching conditions are = 0.1181, Mz = 91.1882 GeV; 

(ml), mb = 4.3 GeV; me = 1.3 GeV. 
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Figure 6: The dependencies of Aan , on the matching point f^b for the 1 — 4-loop 

order cases. For the analytic coupling the curves go above the corresponding curves for 
the perturbative coupling. The normalization and matching conditions are = 

0.1181, Mz = 91.1882 GeV; 
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Figure 7; The dependencies of aan{M^), on the matching point fj,b for the 1 — 

4-loop order cases. For the analytic coupling the curves go lower than the corresponding 
curves for the perturbative coupling. The normalization and matching conditions are 
afo/=5)(M|) = 0.1181, Mz = 91.1882 GeV; 












Figure 8: The analytic and perturbative couplings aaniQ"^), for the 1 — d-looj 

order cases. The normalization conditions are = 0.1181, Mz = 91.1882 GeV 

rrih = 4.3 GeV; , rric = 1.3 GeV. 





